Dipolar interactions in magnetic nanowires aggregates. 
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We investigate the role of dipolar interactions on the magnetic properties of nanowires aggregates. 
Micromagnetic simulations show that dipolar interactions between wires are not detrimental to the 
high coercivity properties of magnetic nanowires composites even in very dense aggregates. This 
is confirmed by experimental magnetization measurements and Henkel plots which show that the 
dipolar interactions are small. Indeed, we show that misalignment of the nanowires in aggregates 
leads to a coercivity reduction of only 30%. Direct dipolar interactions between nanowires, even as 
close as 2nm, have small effects (maximum coercivity reduction of ~15%) and are very sensitive to 
the detailed geometrical arrangement of wires. These results strenghten the potential of magnetic 
composite materials based on elongated single domain particles for the fabrication of permanent 
magnetic materials. 
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I. INTRODUCTION 

Taking advantage of the shape anisotropy for the fab- 
rication of high coercivity materials has been considered 
in the past. We can mention the development of elon- 
gated single domain (ESD) particles PUS] or the AINiCo 
materials [3] in which elongated magnetic structures are 
formed by a proper metallurgical process. These materi- 
als are however limited in their performances to coercive 
fields of the order of a few mT. During the last decade 
it has been shown that elongated magnetic nano-objects 
with well defined shapes could be synthesized jSHS]- We 
have considered the use of magnetic nanorods for the fab- 
rication of high magnetic energy polymer composite ma- 
terials. We have shown that rather high coercive fields, 
up to 600kA/m (= 0.75 T) could be obtained [ID]. How- 
ever, the nanowires are not optimally dispersed in the 
matrix material and tend to form aggregates (Figure 1). 
The scope of this communication is to address the role 
of the dipolar interactions in nanowire aggregates on the 
coercive field in our composite materials. 

The role of dipolar interactions in assemblies of mag- 
netic nanoparticles has already been addressed in several 
situations. In some geometries, such as planar arrays of 



magnetic nanowires fabricated in porous alumina mem- 
branes, the dipolar fields play a very important role. It 
has been experimentally demonstrated that the demag- 
netizing field Hd is simply proportionnal to the packing 
density P: H d = -P.M S [HJ [IS]. Since P can reach 
values up to 30%, the demagnetizing fields are very large 
and significantly reduce the coercivity and more generally 
the anisotropy properties of the arrays [121 [T5] . 




Figure 1. Scanning Electron Microscope images of Co 
nanowires deposited under a magnetic field on a Si substrate, 
(a) At the microscopic scale, the objects form microstructures 
with cigar shapes, (b) Inside these microscopic structures, the 
Co nanowires are rather well aligned with an angular disper- 
sion of the order of a = 7°. 
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The effects of dipolar fields have also been considered 
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for granular hard magnetic materials [14 . In this case, 
Henkel plots |15| -obtained from isothermal remanent 
magnetization and dc demagnetization experiments- al- 
low to put into evidence interactions between magnetic 
grains [16 . In the above situations macroscopic demag- 
netizing field effects must be taken into account. In our 
composite materials, the situation is qualitatively differ- 
ent since elongated cigar shape aggegrates are formed. 
Thus macroscopic demagnetizing field effects, within a 
cigar, are expected to be negligible. In this communi- 
cation we focus on local dipolar field effects between in- 
dividual wires and we use micromagnetic simulations to 
quantify the role of the magnetic interactions between 
nanowires on the coercive field. 



II. MAGNETIC NANOWIRES AND 
AGGREGATES 

We consider the experimental case of Co and CoggiV^o 
nanowires synthesized via the polyol process [5HH] ■ These 
nanowires exhibit diameters in the range 7 — 15 nm and 
lengths in the range 100—250 nm. In the previous studies, 
these nanowires, either randomly oriented or well-aligned 
in solvents, have been magnetically characterized. The 
experimental hysteresis cycles have been discussed within 
the Stoner-Wohlfarth model [TU]. The reversal mecha- 
nism is usually not a coherent one and it is initiated at 
the wires tips. Micromagnetic simulations put into evi- 
dence the role of the detailed shape of these objects in 
the magnetization reversal [T7]. The chemical synthesis 
of magnetic wires has been optimized to achieve a cylin- 
drical shape so as to optimize the shape anisotropy [7] . 

When the nanowires are dispersed in solvents and dried 
under a magnetic field, they tend to form microscopic 
aggregates (Fig. 1) so that there can be direct magnetic 
dipolar interactions between the objects. The nanowires 
are covered with an oxide shell which is typically 1 — 
2 nm thick [TH] so that the minimal distance between the 
ferromagnetic cores is of the order of 2 — 4 nm. 

Figure 2 represents the dipolar stray field at the tip of 
a nanowire (L = 100 nm, r = 5nm, M s — IT) calculated 
using FEMM [T5]. It illustrates that the dipolar field is 
very localized around the tip (in a volume with a typical 
size given by the radius of the nanowire) . The stray fields 
are of the order of 0.1 T at a distance of 4-7nm from the 
nanowire tip (see Fig. 2b). Thus, since the nanowires are 
separated by distances which can be as small as 2 — 4 nm, 
the dipolar field radiated by one wire on its neighbors can 
still be of the order of a fraction of one Tesla. 

In the present paper, we address the question whether 
global or local stray fields may significantly affect the 
magnetic properties of nanowire aggregates. Several ef- 
fects may play a role in modifying the overall coercivity 
of nanowires assemblies. They can be classified in the 
following way: 

(i) The wires may be disoriented with respect to each 
other. This leads to a significant loss in the macroscopic 
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Figure 2. Dipolar field outside the tip of a [L — 100 nm, r = 
5 nm, M s = 1 T) magnetic nanowire (gray area) calculated 
using FEMM |19| . (a) Axial mapping of the induction outside 
the wire, (b) Magnitude of the induction along the wire axis 
and the radial axis. The dipolar stray fields drop quickly but 
are still of the order of 0.1 T at a distance of 4 — 7nm of the 
wire edges. 



coercivity. 

(ii) The fact that the wires are assembled in cigar 
shaped structures has in principle the effect of decreasing 
the macroscopic demagnetizing field effects and increas- 
ing the overall coercivity. 

(iii) The local stray fields created by a wire tip may 
create a local nucleation point and promote the reversal 
of a neighboring wire. 

(iv) The presence of "defective wires" with low coer- 
civity because of a poorly defined shape, a strong mis- 
orientation, or a structural defect may play the role of a 
nanoscopic nucleation point inside the cigar microstruc- 
ture. 

In the following, we quantify these different effects in 
order to estimate the requirements for obtaining compos- 
ite materials with optimum coercivity properties. 



III. MODELLING 

The model objects which are considered are cylindri- 
cal objects of length L = 100 nm and radius r = 5nm. 
The magnetic parameters used correspond to hep cobalt 
epitaxial thin films [20 , saturation magnetization Ms = 
1400 kA.m -1 , exchange constant A = 1.2 X 10~ n J/m. 
The micromagnetic calculations have been performed us- 
ing the Nmag package [21 . The objects were meshed 
using Netgen |22| with meshes consisting of about 1000 
points so that the distance between two nodes is of the 
order of the exchange length £ ex — y // 2A/ /iQMg 3.2nm. 
Increasing the finesse of the meshes did not change the 
results. Magnetocrystalline anisotropy was not included 
in the calculation on purpose so as to avoid masking the 
dipolar effects by extra anisotropy effects. Adding an ex- 
tra magnetic anisotropy to the wires would not change 
the conclusions. 
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A. Orientation effects 

For perfectly aligned infinite ellipsoids, the expected 
coercive field (obtained from the shape anisotropy) is 
M s /2 that is 700kA/m (= 0.88T) in the considered ex- 
perimental situation. However, in the case of nanowires 
of finite length instead of infinite ellipsoids, a significant 
drop in the shape anisotropy is expected (~25%)[T7]. On 
top of that, the misorientation of the nanowires induces 
a loss in the effective coercivity. We assume a Gaussian 
distribution. Let a be the width of the distribution of 
the wires orientations with respect to the applied field. 
H c drops from 0.38 Ms for a very narrow orientation dis- 
tribution (it = 1°) down to 0.21 Ms for a random orienta- 
tion distribution (see Fig. 3). In the observed situation 
of Figure 1, the angular distribution of wires has a width 
of about (7 = 7° leading to a moderate loss of coerciv- 
ity (~ 15%) compared to a = 1°. Thus extremely well 
aligned nanowires do not significantly boost the perfor- 
mances of the material. 



have a large aspect ratio so that we are very close to the 
case of infinite cylinders. The demagnetizing coefficient 
N of a ellipsoid of aspect ratio 10 is equal to 0.02 which 
is very close to the value N = for an infinite cylin- 
der. Increasing further the overall aspect ratio does not 
induce measurable effects. When the magnetic field is 
applied at an angle (ranging from 0° to 20°) with respect 
to the chains, no measurable change in the coercivity are 
observed between the Al, A2 and A3 configurations. 
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Figure 4. Wires configurations used in the calculations. (Al) 
A single wire which is the reference situation. (A2) A chain 
of 2 wires separated by 2 nm. (A3) A chain of 3 wires. 
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Figure 3. Hysteresis cycles for different wires orientation dis- 
tributions: a = 1°, 10°, 20°, 40°, oo(random). In the experi- 
mental situation of Figure 1, a has been estimated at 7°. 



B. Demagnetizing field effects 

The assembly of the wires into cigar shape structures 
potentially modifies the overall demagnetization field dis- 
tribution. To account for this effect we have compared 
the following 3 situations: (Al) an isolated wire, (A2) a 
chain of 2 wires separated by I = 2 nm, (A3) a chain of 3 
wires (Fig. 4). The separation of 2 nm which is a min- 
imum value, corresponds to the experimental situation 
because the wires are surrounded by a shell of non mag- 
netic material which prevents a direct contact between 
the wires. The reference value of the coercive field for 
an isolated wire is He ~ 375kA/m. The micromagnetic 
simulations show that no measurable change in the co- 
ercive field is observed when the wires are arranged in 
chains. This could be expected since the wires already 



C. Effect of local stray fields 

We have considered several situations in which the 
stray field of a wire tip may alter the reversal of a neigh- 
boring wire. We have considered 2 main situations: (i) 
the case where the stray field acts at a wire tip (Bl - 
B2), (ii) the case were the stray field acts on a neighbor- 
ing wire edge (B3). The different configurations are illus- 
trated on Figure 5. We first focused on the situation of 
aligned nanowires before considering misalignment which 
is recurrent in aggregates. The angle between the two 
nanowires long axes is noted a. 

Figures 6(a), (b) and (c) respectively display the hys- 
teresis cycles of the Bl, B2 and B3 configurations for 
different values of the a angle. For the Bl and B2 con- 
figurations with a — 0, the hysteresis cycles are identical 
to the one of a single nanowire and He ~ 375kA/m (Fig. 
6a and b). It implies that for these configurations -close 
to A2- there is no measurable effect of the local stray field 
on the wire tips. Both wires switch at the same time. The 
situation differs for the B3 configuration where the stray 
field acts on a neighboring wire edge. In this case, the 
coercivity is lowered to He ~ 325kA/m, that is a loss 
of 13%, while the cycle remains almost perfectly square 
(Fig. 6c) . The magnetization reversal of both nanowires 
is simultaneous. 

We now investigate the effects of misorientations (a ^ 
0) for the configurations Bl, B2 and B3 in order to repro- 
duce the situation of aggregates. In the Bl configuration, 
the coercive field is reduced from He ~ 375kA/m for 
a = 0° down to He ~ 255 kA/m for a = 40°. The magne- 
tization reversal of both nanowires remains simultaneous 
demonstrating that there is a strong interaction between 
the wires. However in the limiting case where a = 90° 
{i.e. the nanowires are perpendicular to each other), the 
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Figure 5. Wires configurations used in the calculations. (Bl) 
Two wires touching at their tips but shifted by R. The relative 
orientation of the wires can vary from a — 0° to 40°. (B2) 
Two wires touching at their tips but shifted by 2R. (B3) A 
wire creating a stray field around the center of another wire. 



coercive field is somewhat recovered to He ~ 330kA/m 
but the two wires reversals are disconnected. The mag- 
netization of the wire at 90° simply follows a linear de- 
pendance with respect to the applied field. 

The situation slightly differs in the B2 configurations 
since the magnetization reversal of both nanowires is no 
more simultaneous for a = 20° and 40° showing that the 
interactions between wire tips is very sensitive to geomet- 
rical details. Furthermore, in the limiting case a = 90°, 
the coercivity almost recovers the value of the isolated 
nanowire case, suggesting that there are no more mea- 
surable interactions between the wires. 

The situation is qualitatively different in the B3 config- 
urations. Figure 6(c) shows that for a = 10°, 20° or 40°, 
the magnetization reversal of both nanowires is no more 
simultaneous and the coercivity of the aligned wires re- 
covers a value of He ~ 375kA/m corresponding to the 
case of an isolated wire (Al). Surprisingly, it shows that 
neighboring wires in the configurations B3 have measur- 
able effects on the coercive field only when they are par- 
allel to their neighbor and that as soon as there is a small 
misalignemcnt, there are no more interactions. 

Finally, on the one hand, for perfectly aligned wires 
(a = 0) , stray fields leads to a moderate loss of coer- 
civity (13%) only in the B3 configuration. On the other 
hand, misalignment of the nanowires (a / 0) , induces 
a significant loss of coercivity in the Bl and B2 configu- 
rations while stray fields effects become negligible in the 
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Figure 6. Hysteresis cycles of the (a) Bl configuration for a = 
0°, 10°, 20°, 40° and 90°, (b) Hysteresis cycles of the B2 con- 
figuration for a — 0°, 10°, 20°, 40° and 90°, (c) Hysteresis cy- 
cles of the B3 configuration for a = 0°, 10°, 20°, 40° and 90°. 



B3 configurations. The dipolar interactions on the coer- 
civity are very sensitive to the geometrical arrangement 
of wires. 



D. Effect of nanoscopic nucleation or pinning 
points 

We now consider the case of nucleation or pinning 
points created inside a cigar shaped micro-structure by 
wires which have respectively either a lower coercivity 
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or a higher coercivity. In the case of a wire with a ge- 
ometrical defect, corked or with 3 endings for example, 
its coercivity will be decreased. On the other hand, one 
may imagine that some wires contain structural defects 
which pin their magnetization. For the micro-magnetic 
modelling, we have thus pinned the magnetization Mq 
of one of the wires of the structure in a direction along 
the nanowire long axis (Fig. 7 top). Starting from a 
negative saturation field —H sat , when arriving at rema- 
nence, none of the wires has flipped except for the pinned 
one. This pinned wire can thus be considered as a nu- 
cleation point which will promote the reversal of neigh- 
boring wires when the field is further increased. A loss 
of coercivity is expected. On the other hand, starting 
from a positive saturation field +H satl when arriving at 
rcmanence, the magnetization of all the wires are still 
parallel. When the field is further decreased, the pinned 
wire will not flip and will thus act as a pinning point. 
A slight increase in coercivity is expected. Thus, on the 
calculated hysteresis cycle, the effect of pinning can be 
observed on the left side while on the right side, the effect 
of nucleation is observed. 



isolated wire: the pinned wire has no effect. On the 
other hand, when sweeping the field to positive values, 
the pinned wire acts as a strong nucleation point and 
induces a 40% drop of the coercive field with respect to 
an isolated wire. The configuration C2 consists of 9 wires 
among which one is pinned. The reference situation (9 
free wires) is plotted on Fig. 7 as the solid grey line. The 
coercivity is 317 kA/m. In the C2 configuration, when 
the field is swept to negative values, the coercive field of 
the wires chain is unchanged. There is only a small step 
which appears in the hysteresis cycle which corresponds 
to a stiffer reversal for the Nr. 1 wire (see small sketch 
of Figure 7) . We suggest that the different behaviour of 
wires 1 and 2 is due to demagnetizing field effects. When 
sweeping the field to positive values, the pinned wire acts 
as a strong nucleation point and promotes a first reversal 
of wires 1 and 2 at a field 40% lower than the coercive 
field. The other wires behavior is barely modified. The 
coercivity is slightly increased which can be accounted 
for by dipolar interactions with the flipped row. Thus, 
in conclusion, a pinned wire only affects the behavior of 
its neighbor along a chain but not its side neighbors. 



IV. EXPERIMENTAL RESULTS 
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Figure 7. Various configurations of aggregated wires in which 
the magnetization of one of the wire is pinned. (CI) Chain of 
3 wires. (C2) Staggered rows of wires. 



We performed micromagnetic simulations on various 
configurations of aggregated wires. Let us discuss the two 
situations illustrated in the geometries of Fig. 7 (top). 
In the first case, we look at the behavior of a row of 
3 wires in which the middle one has its magnetization 
pinned. The hysteresis cycles of Fig. 7 show the evolution 
of the magnetization of the free wires. That is we plot 
the average magnetization (normalized to one) of all the 
wires except the pinned one. The reference situation of 
an isolated wire is represented as the black dashed line. 
When the field is swept to negative values, the coercive 
field of the wires chain is unchanged with respect to the 



In order to probe the interactions in granular mag- 
netic materials, Henkel [TS] proposed to compare the 
Isothermal Remanent Magnetization (IRM) and Direct 
Current Demagnetization (DCD) curves. The difference 
Sm(H) = (2.IRM(H) - DCD(H))/IRM(oo) - 1 being 
either zero in the case of non interacting particles |23j . 
positive for interactions promoting the magnetized state 
and negative for interactions assisting magnetization re- 
versal. Note that we are measuring the DCD with an 
initial negative saturation field which is opposite to the 
usual convention, hence the change of sign in the previ- 
ous formula for the DCD contribution. We think that 
this way of measuring the DCD is more natural though. 

We performed IRM and DCD measurements on aligned 
Co nanowires deposited on a Si substrate (see Fig. 1). 
The measurements were performed both parallel and per- 
pendicular to the easy axis of the nanowires aggregates 
(Fig. 8). One can notice that 8m(H) is mostly zero ex- 
cept in 2 small regions during the reversal process. The 
5rn(H) value is small (max ~ -14%) and negative sug- 
gesting that only dipolar interactions assisting reversal 
take place. This is in agreement with the above calcu- 
lations by showing that (i) forming aggregates does not 
stabilize the magnetic structures and that only extra nu- 
cleations can be induced by dipolar effects and (ii) that 
these effects remain rather small even in very compact 
structures with densities close to 1. No qualitative differ- 
ence is observed in the measurements along and perpen- 
dicular to the easy axis. The facts that two small peaks 
are observed in Sm(H) during the reversal process sug- 
gests that different mechanisms are taking place at the 
beginning and at the end of the aggregate reversal. This 
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needs to be investigated further. 
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Figure 8. IRM, DCD and Henkel plots measured on aligned 
Co wires in a polymer matrix. The magnetic field was ap- 
plied either parallel (top) or perpendicular (bottom) to the 
aggregates axis. 



V. CONCLUSION AND PERSPECTIVES 

We have investigated the role of magnetic dipolar inter- 
actions between nanowires in aggregates via micromag- 
nctic simulations. In the case of aggregates composed of 
perfectly aligned nanowires, these simulations show that 
no major change should be observed in the shape of the 
hysteresis cycles or in the values of the coercivity and the 
remanence. In the worst situation a drop of 20% in the 
coercive field might be observed. We underline that local 
dipolar interactions effects on the coercivity are very sen- 
sitive to the detailled geometrical arrangement of wires. 
As for the effect of nanowire misalignment inside aggre- 
gates, we show that even very poor alignment leads only 
to a 30% reduction in the coercive field. Using Henkel 
plots we have experimentally confirmed that dipolar ef- 
fects are indeed small in these aggregates of nanowires or 
that at least they do not significantly modify the mag- 
netic behavior. This validates the way of interpreting the 
experimental magnetic hysteresis cycles in |10| . 

It is striking that in dense aggregates of magnetic 
nanowires, dipolar interactions play a very limited role. 
We conclude that in dense materials based on magnetic 
nanowires, the interactions between wires are negligible 
and that the optimization of the materials goes through 
the optimization of the properties of individual wires. 
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